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ABSTRACT 


Atmospheric turbulence can presently be measured only on the ground 
or by means of airplanes flying through the turbulent area. A properly 
designed radar system may be able to detect and analyze atmospheric 
turbulence from a position remote from the turbulent area. The choices 
of antenna, wavelength, transmitter power, receiver sensitivity, and 
detection methods are examined. The differences between pulse or cw 
and bistatic or monostatic systems are discussed. A comparison with 
existing radar systems used for atmospheric turbulence research is made. 
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I. INTRODUCTION 


The present study was instigated through a grant from the Electronics 
Research Center of the National Aeronautics and Space Administration for 
a study of clear air turbulence and methods for its detection. The present 
treatment includes detection methods of regions of disturbed refractive 
index in a broad sense and is not limited to clear air turbulence alone. 

Radar return from the atmosphere is due to many factors including 
particulate matter such as rain, snow, sleet, bugs, birds, dust, etc. , but, 
in addition, a radar signal may be scattered from regions of disturbed 
refractive index of the air. The scattering cross section of regions of 
disturbed refractive index is much less than the radar cross section of 
particulate matter such as rain. 

The object of this study is to determine the radar system parameters 
necessary to detect and analyze iiadar returns from regions of disturbed 
refractive index, in particular, those related to clear air turbulence. 

The major factors to be investigated are the system wavelength 
dependence, the difference between monostatic and bistatic radars, the 
desirability of using cw (continuous wave) or pulse radar and find the 
overall system performance required to detect and analyze turbulent 
regions . 

A relatively simple model of the turbulent region will be used but the 
most critical features of size and degree of turbulence of the region will 
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be retained. An attempt will be made to state the assumed model clearly 
with references included to support the choice of model. 

A. Nature of Clear Air Turbulence . 

Clear air turbulence is defined as relatively severe turbulence 
of the atmosphere at high altitudes, usually above 20, 000 feet. At lower 
altitudes, the atmosphere has a larger variation of refractive index due 
to the presence of water vapor. This greater variation produces a larger 
radar cross section for the low altitude case. 

The scattering cross section of region of turbulent air will be 
evaluated both for backscatter and for scattering at other angles. In 
practice this is the weakest link in the chain of parameters involved in 
the detection system since our knowledge of the scattering cross section 
of the turbulent atmosphere is very limited. A few experimental results 
can be cited to obtain an estimate of the error involved in the application 
of the theory. 

The radar equation is examined to determine the important factors 
in obtaining a usable return from the region of disturbed refractive index. 
These factors are then discussed individually. 

The most important factor is the general structure of the turbu- 
lent region, which will determine the scattering cross section, and the 


size of the antenna. 
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Two difference types of emission are discussed, pulse and cw, and 
the methods of detection for each are examined. 

B. Uncertainties in Refractive Index Structure and Suggestions 
for Clarification. 

The characteristics of the variations in the index of refraction in 
the atmosphere are not completely understood, especially predictability. 
There have been several papers written on the possibility of detecting clear 
air turbulence by radar techniques. Some confirmation of 

turbulence theory, as applied to scattering of radio waves, is needed. 

Variations in the dielectric constant of the air may or may not 
be easily related to the type of turbulence encountered by aircraft. 
Especially it would be desirable to discover what are the type of variations 
and the structure of these variations for the case of clear air turbulence. 
One way of more closely approaching this objective would be to make 
simultaneous radar measurements and refractive index measurements. 



II. RADAR EQUATION 


The relationship between the power transmitted and the power received 
is known as the radar equation 9 . An examination of this equation will be 

necessary to determine which parameters can be optimized for detection of 
regions of disturbed refractive index of the atmosphere. 

The radar equation is 


p R = 


G_G_ X a VP^ 
T R v T 

3 2 2 

2-(4tt) R r Z R t 


where 

G^ = gain of transmitting antenna 

G = gain of receiving antenna 
R 

X = operating wavelength 

a = scattering cross section per unit volume 
V = effective scattering volume 
P = power transmitted 
P = power received 
= range to receiver 
R = range to transmitter 
Cj = beam filling factor. 

The distribution of power across the beam is not uniform as assumed 
in the radar equation. A correction factor to reduce the power received 
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may be added to more nearly estimate the true received power. This factor 
is derived for a gaussian distribution of power across the beam by Probert- 

# 1 r \ /OQ\ 

Jones' ; and Stephens' . This factor is 2 tn 2 = 1. 39. In order to have a 
more conservative result and because the distribution is not gaussian, a factor 
of 2 is chosen. The various parameters in the radar equation are not inde- 
pendent of each other. For example, the antenna gain and beamwidth are 
related to wavelength. The scattering volume, which is defined by the pulse 
length and antenna beamwidth for the monostatic system or by the intersection 
of the two antenna patterns for the bistatic system, is dependent also on wave- 
length as well as other factors. 

The size of the antenna beam defines the volume of turbulent air that 
will be illuminated. As the gain of the antenna increases, the effective 
scattering volume decreases. Hence, the optimum antenna size and gain 
need to be investigated with respect to changing wavelength. For best 
efficiency, the antenna beam should be filled, requiring that the extent of 
the turbulent region be estimated and the effective scattering volume shaped 
accordingly. 

The standard radar parameters including transmitter power and 
modulation, receiver noise figure and detection systems, will be evaluated 
and conclusions drawn as to the optimum overall system for studying the 
atmosphere using radar. 



III. ATMOSPHERIC TURBULENCE 


Random variations of refractive index of the atmosphere occurring in a 
distance of the order of a wavelength will cause scattering of electromagnetic 
energy. This scattering is very small and thus a very sensitive radar system 
is necessary to detect it. 

A. Turbulence Structure 

The variations in refractive index with distance can be characterized 

( 25 ) 

as a random process with a space spectrum as shown in Figure 1. Tartar ski v 7 
has shown that radar return is a function of the structure function defined by 

D n (r) = [N(o) - N(r)] 2 

where N(o) and N(r) are the refractive indices at two points separated by a 

( 21 ) 

distance r. When r is in the inertial subrange 

2 2/3 

D (r) = C r L « r « L 

n n m o 

2 

where is a constant depending on the strength of the turbulence. For 

r very large, the two samples N(o) and N(r) become independent and thus 

[ N(o) - N(r)] 2 = 2 (An) 2 . 

Letting r equal L^ gives 


, ,2 1 2 _ 2/3 

(An) C n L o ‘ 


A more exact formula can be derived taking into account the fact 

2/3 

that the r variation does not extend to L . This formula is 

o 
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Refractive Index Spatial 
Spectrum 


Figure 1 


8 


( 1 n, 2 ,.19L 2 / 3 C 2 . 

o n 


The one dimensional spatial spectrum that corresponds to this 


structure function is^^ 

F' (k) = . 124 C 2 k' 5/3 . 
n n 

The Obukov spectrum is a normalized version of this. It is 

_ . _ 2 2/3 , -5/3 

F (k) = — k k 
n 3 o 

The radar return is controlled by the amplitude of the turbulence 
spectrum over a narrow range of the spacial wavelength centered at the 
radar's wavelength. This is represented in Figure 1 by the narrow band 
filter. Spectra have been measured over a limited range in the inertial 
subrange Austin, Texas using a set of four UHF (400 mHz) refractometer s 
mounted in a three-dimensional array. With spacings from 30 cm to 4 
meters, the mean square value of difference in refractive index between 
temperatures was measured at a height above ground of about 75 meters 
and found to correpond to the spectrum in Figure 1. The data are, of 
course, nonstationary and contain large bursts probably caused by patches 
of higher water vapor content of the air at that instant. This would make 
the refractive index fluctuations much larger at low altitudes, enabling 
radar data to be obtained with a relatively simple system. 
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B. Clear Air Turbulence 

Turbulence of the type that occurs in clear air lies in layers 500 
to 3000 feet thick and up to 20 miles in horizontal extent. The system used 
to detect this turbulence must take this fact into account. An ordinary 
airborne radar may have only a very small percentage of its beam filled 
at a range sufficient to permit an aircraft time to take evasive action, 

C. Problems in Detection 

A jet plane must have at least one minute's warning in order to 
avoid turbulent areas. For a plane flying 600 mph, this requires a 
minimum detection range of ten miles. Some method of separation of 
the ground clutter from the air return must be employed. Without this 
separation, the radar window (area on the scope clear of ground clutter) 
for a pulse system is determined by the height of the plane which will be 
about six miles for the normal cruising altitudes of subsonic jet planes. 
Without the elimination of ground clutter, an adequate warning will not 
be available to the crew. Another problem with airborne radar is that 
of determining the height of a disturbance. At a range of ten miles, a 
five degree beamwidth (probably the minimum angle for the relatively long 
wavelength which appears to be better for turbulence detection) receives 
radar return from a height range of 5000 feet. This uncertainty in height 
gives the pilot little chance of avoiding the turbulence. 
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Thus, a ground-based research system seems to be the best 
possibility at present. Either a pulse monostatic or pulse or cw bistatic 
system could be used. A system of each type is described in this report. 

As indicated earlier, clear air turbulence of the type that air- 
craft encounter seems to lie in layers 500 to 3000 feet thick and from 1 to 

(8 24 ) 

20 miles in horizontal extent' * . A vertically pointing radar would, 

therefore, have 100 per cent of its sampling volume filled with turbulence 
at some time as the turbulent area passes over. A ground-based bi- 
static system would also work well for the same reason. 

D. Available Data 

The general type of refractive index disturbance in the atmosphere 
up to 75 meters altitude has recently been measured in Austin, Texas with 
a three dimensional array of UHF r efractom eter s . These measurements 
show that the turbulence spectrum is essentially the same as the theoretical 
one with values of about . 25 N units occurring frequently, N = (n - 1) 10^ 
where n is the index of refraction. 

Average values of of 10 ^ cm occur frequently. The 

structure of the turbulence appears to be mainly horizontal eddies. The 
size limit has not been reached with a refractometer spacing of 27 feet which 
was thq maximum obtainable at the time. A reasonable assumption is that 
the eddies continually increase in intensity with size, with reasonably 
homogeneous turbulence at the smaller scales (one centimeter to one meter) 


of interest for radar measurements. 
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Structure function measurements indicate that the turbulent eddies 

2 2 

are substantially horizontal with horizontally greater than C vertically 

by about 25 per cent. Measured C ^ ranges from 10 ^ to 10~^ with a five 

- 14 

minute average of 10 at 80 meters height. At this altitude, water vapor 

has a strong effect on the radio refractive index of air and this factor 

2 

undoubtedly causes the high values of C that were measured. The 

n 

2 

variation of C n with altitude is essentially unknown. 

Variations of refractive index with height have been measured. 

^ 2 

A AN of . 2 may occur up to 16, 000 feet altitude. The AN variation with 
height is not constant but, in general, decreases with increasing height 

(9) 

except for turbulent layers. 

The type of instruments used in obtaining this data is described 

in reference 18, the first section. 

(24) 2 

Stephens and Reiter' 9 estimate the magnitude of C for moderate 

n 

turbulence. Atlas, Hardy and Naito^ estimate C ^ to be 10 ^ cm 

n 

-14 -2/3 (2) 

for moderate CAT and 10 cm for the most reflective CAT. ' 1 For 

the case under consideration, that is, a ground based radar observing 

relatively low altitude (up to 15, 000 feet) turbulence, the magnitude of 
2 

C^ should be larger. Also, it is thought that water vapor may play a more 

significant part in the measurements made here which might account for the 

2 -13 2 -15 

larger C^ measured here (10 > C^ > 10 ). The purpose of corre- 


lating refractive index measurements with radar measurements to investigate 
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the general properties of atmospheric turbulence will be served by low 

altitude measurements and thus will not require as powerful a system 

as expected for Clear Air Turbulence. 

Buehler and Lunden of the Boeing Company have had fair success 

with a VHF system (200 mHz) in detecting turbulence^’ ^ * 2 ^, They 

have installed a system in a Boeing 727 and flown it for 60 hours and 

have nine confirmed turbulence encounters of both radar sighting and 

airplane vibration at high altitude in clear air. 

The radar volume scattering cross stection for atmospheric 

(2) 

turbulence has been derived from Tatarski's work by Atlas, et al v , 

(24) 

and Stephens and Reiter for random variations; for periodic variations 

( 22 ) 

by Smith and Rogers' ; and for a uniform gradient by Buehler and 
Lunden^ . 

E. Scattering Cross Section 

The most useful scattering cross section in this case is the 
one based on Tatarski's theory. The general scattering cross section 
per unit volume is 

tt sin 2 p (An) 2 R 2 F (K) 

n 

a = - 4 e 

8 sin — 


where 


0 = angle between the transmitter and receiver beams 
(scattering angle) 
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(3 = angle between the direction of the receiver and the electric 

field of the transmitted wave 

2 

An = mean square refractive index value 

2 

F (K) = one dimensional spectral density of An 
n 

4tt . 0 

K - T sin- 


2 2/3 -5/3 

The Obukhov spectrum is F (k) =*? k k where k < k < k . 

n 3 o o m 

k = - 7 “ and k - 7 -- — where L and L are the inner and outer scales, 
o L m L mo 

o m 

respectively. L is approximately 1 cm for moderate turbulence and L 
m o 


is from 100 to 600 meters 


(2, 19) 


The F^(k) for the calculation of the scattering cross section is 

along the dotted line in Figure 2. The scattering cross section, in terms 
2 

of C , is derived below, 
n 


tt sin 2 g (An ) 2 Y? F (K) 

n 


a = 
v 


Q . 4 e 

8 sin - 


_ , TT sin 2 0{&n) 2 K 2 , ,2 „ 2/3 „-5/3, 

°v' ( „ '.Ve )( 3 K o K ’ 

8 sin ^ 


2 2/3 2 

(An) =0.19 L C 
■ o n 


.2/3 


,1/3 


. 2 


TT x 2 /3 _ 2 n r (2 tt) r (4rr) ' . 9 .1/3, sin 

= — [ 0 . 19L^ ] V s m — ) ] 


12 


o n-- L 2/3^ x l/3 
o 


2 ' J . 4 
sm 


col CD | CM 
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covers most of the usable radar spectrum with scattering cross section 

peaking at X = 5 or 5 cm. ^ The X ^ law is only accurate for 

X > 10 L . 

m 




IV. ANTENNA CONSIDERATIONS 


A. Volume for the Monostatic Case 

The scattering volume will be calculated assuming the radiated 
power is concentrated within the 3 db points of the antenna beam. There 
is an additional effect which must be taken into account for backscatter 
of a pulse from a continuum of particles. At any instant after steady 
state, the power returned is due to the number of particles occupying 
the volume due to one half the pulse length (Figure 4). 

This effect is due to the two way travel of the radar pulse. If 
the scattering angle 0 is less than 90 degrees, the factor of one half can 
be removed. 

The scattering volume must now be calculated (Figure 5). 


i ^ 


2 " I 3 R+ I 


= J J J r sin 0 dr <36 dcp = J J [ ^ sin 0 d8 dcp 


o o R-- 


o o 


R -2 


3 3 

2tt r . h. . . h. -i r . of 

- 3 - [ ( R+- J ) - ( r - 2 * ] [ 1 - cos - ] 


„ , 3 ^ ,3 2 4 

2 n r h r h r i /i ^ , Qf 

^ L ^ 4 3R h] [ 1 - cos '^‘J — ^ L ^ "1 3R hj ^ "^4 


.)] 


IT r h ot-,2. 2 

= - [— + 3R h] a . 

O' = . 612 — radians where a is the diameter of the aperture. 
a 


1 ? 
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With this substitution, 

h 3 2 > 3 

V = . 392 [ — + 3R h] ( - ) . 

4 J a 

The pulse length will be chosen to optimize the returned signal 
from a horizontal layer 500 to 3000 feet thick. Therefore, for a 1000 
foot thick layer (a median value), aptil.se length of 2ps will be optimum 
as h = * 2 “ 

B. Volume for Bistatic Case 

For the bistatic radar case, the scattering volume will be 
assumed to be defined entirely by the intersection of the antenna beams. 
The scattering volume geometry is illustrated in Figure 6. The 
problem of calculating the volume can be much simplified if the beams 
are assumed cylindrical. This assumption will hold when the length 
of the intersection is less than of the range to the volume. The two 
intersecting cylinders are shown in Figure 7. 

The scattering volume may be thought of as a stack of parallel- 
ogram shaped slices of dz thickness. The sum of all these slices equals 
the total volume. 

To calculate the common volume, take slices of the intersection 
parallel to the xy plane (Figure 8). Next, take a cross section parallel 
to the yz plane to get the variation in c and d which are the edges of the 
parallelogram (Figure 9). 


GROUND 


Scattering Volume for Bistatic Radar 


Figure 6 




Model for Calculating Scattering Volume 
for Bistatic Radar 


Figure 7 
Y 



Area of parallelogram = cd sin cp 
Par allelgram 


Figure 8 
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2 2 

z ,y .2 

— + — s,n = 1 

a a 


is the equation of this ellipse (Figure 8). From this 


y =- 


2 

V a - z 
sin 


Now take cross section parallel to the xz plane (Figure 10). The 

equation of this circle is 

2 , 2 2 

x + z - b . 


From this 


= vh 


,,2 2 

x = V b - z 


Sum up the volumes of all the parallelograms of depth dz, dimensions 

Jz 2 _rz 2 

. 2v a - z , 2V b - z 

d : — and c = 


sin 0 


sin 0 


V = 2 J 

o 


„ /~2 2 r 2 2 

4v a - z vb - z 

“ 2 . 
sin, 0 


sin 0 dz where d = min {a,b}. 


For the case when a = b 


V = — 


8 


r , 2 2 % , _ 16 a' 

(a - z ) dz - - — : 

J 3 sm cp 


16 a" 


sm cp 


o 
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Since cp - 180° - 0, sin cp - sin 0 and the volume becomes 
16a 3 

v = 3liT? ‘ 

But a = R-^- 

Cd 


where at is the antenna beam angle and R is the range to the center of the 
volume. So the volume is 


V = 


16 (R f> 


3 3 
2 R a 


3 sin 0 3 sin 0 


= . 153 


3 3 

— ( A ) 

sin 0 v a 


in terms of parameters already defined and used in the radar equation. 

C. Limitations on Antenna Size 

The size of the antenna determines both the gain and the scattering 
volume, so that the choice of antenna size is not immediately obvious. 
Another factor of importance is the change in the antenna characteristics 
at close range due to the transition from the Fraunhofer to Fresnel regions. 
This gives a minimum range limitation which is 



R = range 
D = diameter 
\ = wavelength. 

For this criterion, the true gain is . 94 of the Fraunhofer gain which is a 

( 20 ) 

small loss for this application . Page 28 gives a more complete analysis 


of this situation. 
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There are two other practical limitations which are the maximum 
physical antenna size and the maximum obtainable gain. The physical size 
of the antenna is purely a matter of choice and cost at long wavelengths but 
at short wavelengths the tolerance required on the surface of the reflector 
makes antennas with gains greater than 60 db impractical for general use. 

A practical limit would be 50 db for an antenna gain that can be obtained 
without too much difficulty. The antenna size will be limited by this factor. 

For longer wavelengths, the antenna size is limited by structural 
considerations. The size and weight of the antenna present prohibitive 
requirements for the supporting structure and positioning mechanism. A 
reasonable choice of antenna size is 15 meters (15 feet) diameter. 

An approximate formula for the gain of a parabolic antenna with 

tapered illumination at the edges, 10 db (power) below the intensity at the 

, . (26) 
center is 



A = aperature area ) 

) same units 
X = wavelength ) 

The cross over wavelength can be found by solving for the wavelength 


at which a 15 meter antenna has 50 db gain. 
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At wavelengths longer than 11 cm, a constant antenna size of 15 
meters will be assumed and at wavelengths shorter than 11 cm, a constant 


gain of 50 db will be assumed for the maximum antenna size and gain. 



V. SYSTEM WAVELENGTH DEPENDENCE 

The system effectiveness as a function of wavelength can now be found. 
The radar equation is 

_ G T P tV CT v V 

R 2 (4tt ) 3 R r 2 R t 2 1 

where P , R D , R will be assumed constant with respect to wavelength for 
T R T 

this calculation. The scattering cross-section is 

, Q , , . -.2 9.- 11/3 -1/3 _ 2 

= . 394 (sin g) (sin-) X . 

-1/3 

o has a X wavelength dependence for wavelengths greater than about 

-1/3 

10 cm but does not obey the X law for shorter wavelengths. For 

( 2 ) 

wavelengths between 10 cm and 3 cm will be approximately constant. 

A. Backscatter Case 

First, the radar equation for the single antenna, pulse radar will be 
written. For a lps pulse length, the volume becomes 


\ 3 2 3 

V = 35 ( — ) [7.5x10 + R ] meters 

a 


R »a 


where a = diameter of antenna. 


For wavelengths shorter than 11 cm, G_, G and a will be constant, 

1 R v 
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p - . 2 V ( P T G T G R CT v °1 

R 2 D 2^ 2 ) 
2(4n) R R 


, , 2 , P_G_G_a C, 

= \ 2 [353 {|) ( 7.5 x 10 3 + R 2 )] ( — - 

2(4tt) R R 


= (constant \ . 

This equation shows that for wavelengths shorter than 11 cm. there is a 
severe penalty in received signal to noise ratio for the assumed system 
as the wavelength decreases. 

However, for wavelengths greater than 11 cm, the situation is 
quite different. The gain is no longer constant and neither is a . The 
gain is 



The received power now becomes 


P T C 1 

P r = < i 2 2 > 

2(4tt) R r R t ^ 


x ] [X 2 ] [• 394\" 1/3 C 2 ] [353 £) (7. 5 x 10 3 + R 2 )] 

. o n a 


= (constant ) X 


1/3 
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The received signal strength is now substantially independent of 
frequency with a slight penalty for longer wavelengths. Over a wavelength 
range of 10 cm to 1 meter there will be a reduction in signal-to-noise ratio 
of only about 3 db which is fairly small. However, good receiver noise figure 
and high transmitted power are more easily obtained with less cost at the 
longer wavelengths. 

In order to have a small sample size to get some idea of the 
statistics of the turbulence and to have a symmetrical volume, the wave- 
lengths between 10 cm and 30 cm would probably be the best. Figure 11 
shows a plot of the relative efficiency of the system against wavelength for 
the backscatter case. 

At the longer wavelengths and long ranges, the volume may be- 
come so large as to not be completely filled causing a drop in the efficiency 
greater than that shown. If the curve for constant antenna size were extended 
for X < 11 cm, it would level off due to the change in the dependence of on 
wavelength. 

B. Bistatic Case 

For the bistatic radar, the wavelength dependence will be different 
due to the fact that the volume is proportional to the third power of the beam 
width instead of the second power. 

For scattering angles (0) less than about 120 degrees, the radar 
wavelength will appear to be lengthened. The radar samples the magnitude 
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of the structure function F (K) at a wavenumber K = sin (ry) (seeFigure 1); 

therefore, as the scattering angle 8 becomes less, the wavenumber K becomes 

-1/3 

less. The effect of this is to cause the X law for the scattering cross 

section to apply for shorter radar wavelengths. 

The wavelength dependent for the bistatic case will be found using 
the same criteria for the antennas of a constant gain of 50 db for X < 1 1 cm 
and a constant diameter of 15 meters for \ > 11 cm. 

For the short wavelength case, the radar equation becomes 


P R = 


P T G T G R C 1 X V<T y 

2(4tt) 3 R 2 R 2 
R T 


/ „ .x „2 V . R 3 .X 3 , ... . 2 . -11/3 9 - 1 / 3_ 2 

= (constant) (X ) (. 153 — : — - (— ) ) (. 394 sin 8 sin — X C ) 

sm y a £ . n 


, 4. ,,,2 3-1/3 

= (constant) XXX 


14/3 

(constant) X , X < 11 cm. 


For the long wavelength case 


P R 


P T G T G R C i x V °v 
2(4tt) 3 R r 2 R t 2 


= (constant) G„G X Va 
1 i\ v 


, „ ,7rra 2 , „ 2 W lc „ R 3 3 V , „„„ . 2„ . -11/36,-1/3^2, 

= (constants — - — ) (X )(. 153 — ■ -■ ■(—) ) (. 394 sm (3 sm — X C ) 

. c, smu a L n 

X 

2/3 


= (constant) X 
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The reason for the greater dependence on wavelength for the bistatic 
case is that the volume is defined completely by the beamwidth whereas in the 
monostatic case the pulse length defines one dimension. These results are 
plotted in Figure 12 for comparison with the backscatter case. Notice that 
the curve for \ > 11 cm rises instead of drops. 

For low altitudes (less than one or two kilometers) disturbed 
refractive index, the scattering volume will be small enough that it will be 
fairly well filled and these system efficiencies will be realized. 

There will be a minimum range for which radar studies can be 
carried out. This is due to being within the Fresnel region of the antenna 
for which the basic assumptions of the radar equation are invalid. For both 
the bistatic and backscatter radars, the range will be that from the scatter- 
ing volume to the antenna but the bistatic system will be able to look at lower 


altitudes and still measure the F^(K) vertically (see Figure 2). 

The boundary between the Fraunhofer and Fresnel regions is 

d 2 

actually a gradual change from one to the other. By using R = -r— as 

mentioned before, the predictions of scattering volume and gain of the 

antennas will still be fairly close to the true values. From Figures 13 and 

d 2 

14 the peak radiated power density occurs at R = . 2 —— . The beamwidth at 


this range is approximately twice that in the far field. A plot of the 


minimum operating range vs. wavelength is given in Figure 15 for the 
previous assumptions; that is, for X < 11 cm, a constant antenna gain of 





HO ID 






Dependence of the Minimum Range on Wavelength 
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50 db is assumed and for X > 11 cm, a constant antenna size of 15 meters is 
assumed. 


For accuracy in gain and beamwidth calculations, the minimum 

d 2 

range — — should be used; however, the minimum range at which usable gain 

X 2 
. 2D 

is obtained is at R =— but inaccuracies in the antenna surface and phase 

error across the aperture make, the gain and beamwidth calculations only 
approximate. 

For low altitude observations with a 50' antenna, the best wave- 
length would be about 30 cm giving a minimum range of about 750 meters. 

At low altitudes, it would be possible to measure with refractometer s the 
disturbed refractive index at the same time as the radar measurements are 
taken. The antenna size would have to be reduced below 15 meters for 
shorter wavelengths in order to stay within the far field when the backs catter 


system is used. 



VI. RECEIVER 


A very sensitive receiver will need to be employed in order to achieve 
the maximum signal- to- noise ratio on the radar return. The ordinary 
crystal mixer front end is far too noisy unless extremely high transmitter 
power is used. This consideration dictates the use of some sort of low noise 
amplifier before the mixer. There are several types available; among these 
are tunnel diode amplifiers, parametric amplifiers using varactor diodes, 
masers using ruby crystals, traveling wave tubes, and transistor amplifiers. 

A. Signal- to-Noise Ratio 

The noise level at the input to the receiver limits the usable re- 
ceiver sensitivity. This noise comes from several sources for microwave 
radar. The most important of these is the noise introduced by the receiver 
itself. For example, the equivalent noise temperature of the sky at micro- 
wave frequencies is approximately 50°K, but the receiver noise temperature 
for a crystal mixer (usual type) will be about 2500 to 5000°K. 

The equivalent noise temperature is defined as the temperature of a 
resistor at the input of a noiseless receiver (matched to the input impedance 
of the receiver) that produces the same noise power output as from the noisy 
receiver. The relationship used is 

N = K T B 
b e 
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where 

N = noise power 

-23 

K, = 1.38 x 10 (Boltzman's constant) 
b 

B = integrated noise bandwidth 

T = effective noise temperature in °K, 
e 

The noise power output is independent of the resistance which makes the 

power relation more useful than the voltage one since 

V XT = 4K T BR 
N be 

and the receiver input impedance must be known. The bandwidth is the 

integrated noise bandwidth which is 
00 

B = J H(f)df 3db bandwidth. 


A more useful measure of the effectiveness of the system than the 
received power alone is the received signal-to-noise ratio. 


S 

N 


received = 



where 


Now, 


N = noise power introduced in the system. 


S 

N 



P 


R 


K, T B 
b e 


(assuming white noise). 
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T is defined the following way: 


(F -1) (F -1) 

T = T + T . [(F.-l )+— £ + —4; + . . . +„ _ 

e sky ref 1 G r G,G, G,G„. . .G 


(Fj-1) 


■ + . . . ] o 


1 


12 


1 2 


i-1 


K 


where the noise figure F. is the signal- to-noise ratio at the output of stage i 
divided by the input signal-to-noise ratio, G^ is the gain of stage i and T ^ 
290°K. The radar equation now becomes 


S . P T G T V V C 1 

N „ , . . 3_ 2_ 2_. _ _ 
2(4tt) R r R t K b T e B 


B. R. F. Amplifiers 

All of these types of R. F. (radio frequency, i. e. , before the mixer) 
amplifiers are not suitable for this application. 

The maser amplifier is very expensive and requires a complicated 
cooling unit. For this application, all of the leakage power cannot be elimi- 
nated. The maser saturates with a power level greater than about -35 dbm 
which is not too difficult to obtain. However, it recovers very slowly (1 ms) 
which makes it almost impossible to use for a pulse system. These objections, 
along with the cost, make the maser a poor choice for this purpose. 

Traveling wave tubes are expensive and do not exhibit better noise 
figures than par ametric amplifiers and tunnel diodes. Therefore, these 
would probably not be as useful. 
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For microwave frequencies (X < 30 cm), this leaves the parametric 

amplifiers and tunnel diode amplifiers. Both of these have fast recovery 

(4) 

time (< 1 M* s ) and both have high saturation levels' . 

Parametric amplifiers can have noise figures of about 2 db while a 
noise figure of 4. 5 db is typical for a tunnel diode amplifier. Either ampli- 
fier will be suitable for use in a turbulence detection system. 

At wavelengths about 1 meter, a greater variety of amplifiers is 
available. Transistor amplifiers can be built for low cost with noise 
figures of 2 db for frequencies of 400 mHz and below. Thus, greater 
receiver sensitivity with less cost is possible at VHF frequencies. A plot 
showing the relative noise figures of various receivers as a function of 
frequency is given in Figure 16. 

C. External Noise 

Another factor limiting the effective noise temperature in this 
frequency range is radio noise from extraterrestrial objects and the earth's 
atmosphere. This noise is a function of frequency and the pointing direction 
of the antenna. Noise from extraterrestrial objects can be minimized by 
pointing away from the galactic center, sun, and moon. The atmospheric 
noise is minimum when the antenna is pointing vertically and is maximum 
for horizontal pointing. 

The maximum and minimum sky temperature as seen by an ideal, 

( 21 ) 

antenna are given in Figure 17. The maximum and minimum curves 
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Receiver Performance 


Figure 16 




SKY TEMPERATURE ( °K ) 
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FREQUENCY (m Hz) 


Sky Temperature 


Figure 17 
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for galactic noise occur when the antenna is pointed at the galactic center 
and the galactic pole. The maximum and minimum curves for atmospheric 
noise are for horizontal pointing and vertical pointing. 

The atmospheric and galactic noise set the minimum usable noise 
figure of the receiver to about 10 to 50 degrees K. 
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VII. SPECIFIC SYSTEMS 


Two types of radar systems for research on disturbed refractive index 
in the atmosphere acetd be designed. The first type is a monostatic pulse 
radar with a vertically pointing antenna. The second is a bistatic radar 
with the two antennas widely enough separated to give a scattering angle 
between 90 and 160 degrees. The system parameters are estimated for 
X = 3, 10, 30, and 100 cm. 

A. Monostatic (backscatter) Radar 
1 . Radar Equation 

The radar equation for this case becomes 


S 

N 


2 2 

G * P T°v 
128rr 3 R 4 K b T 


B 

e 


x 


VC 

2 


where 


— = signal-to-noise ratio at the receiver output 

R = range 

G = antenna gain 


-1/3 2 2 -14 -2/3 

= ^ 394X C , C ranges from 10 cm at low 

n n ° 


K, 


~16 2/3 

altitudes to 10 cm or less at higher altitudes 
= 1. 38 x 10" 23 joule /°K 


B = receiver bandwidth 


h = pulse length 
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= power transmitted 
X = wavelength 

C = beam filling factor 

3 Z 

li Z % “ 

V = . 392 ( — + 3R h) ( — ) = scattering volume 

tc a 

2. System Parameters 

The pulse length is chosen as 2p,s (600 meters in space) so 

that the scattering volume is approximately as long as a turbulence layer 

is thick. The minimum receiver bandwidth is 1 mHz due to the pulse 

length of 2p.s. The beam filling factor will be approximately unity for 

most clear air turbulence as it lies in horizontal layers about 300 meters 

thick and up to 10 km in extent. The receiver RF amplifier has an 

assumed noise figure of 4 db for X = 3 and 10 cm, 3 db for X = 30 cm, 

and 2 db for X = 100 cm. The gain is assumed to be 15 db. Either 

tunnel diode or parametric amplifiers can meet these specifications. 

The amplifiers for X = 30 or 100 cm could be a transistor amplifier which 

would be less expensive than the parametric amplifier. The receiver 

mixer is chosen as a crystal mixer with a noise figure of 9 db. The 

IF amplifier has a 1 me bandwidth with a noise figure of 2 db. Losses 

in the transmission line and duplexer will be assumed to be — db. The 

antenna noise temperature is assumed to be 60°K. This noise is due to 

( 21 ). 

the sky temperature plus noise entering via the side lobes . The 
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system parameters are listed in Table 1. Figure 18 shows the single pulse 
received signal-to-noise ratio as a function of range for the four frequencies. 

Figure 18 shows that the performance of the 3 cm system is 

much less than that of the longer wavelength systems. This is largely due 

to the reduction of the effective area of the antenna for receiving. Figure 

18 also shows that at these short ranges, only a moderate peak power of 
5 

10 watts would be necessary for a reasonable signal-to-noise ratio in 
order to analyze the type of disturbance causing the return. It is desirable 
to limit the sample size in order to be able to tell more about the structure 
of the turbulence which would be best accomplished with either the 10 cm or 
30 cm systems. As far as cost is concerned, the 100 cm system would be 
the optimum. 

3. Detection Methods 

The method of detection can change the effective signal-to-noise 
ratio of the received signal from that predicted by the single pulse signal- 
to-noise ratio to which the radar equation as used so far applies. 

The maximum rate of change of the signal scattered from a 
turbulent region will be about 50 Hz to 100 Hz. This means that each pulse 
will not be independent of the preceding pulse if the prf (pulse repetition 
frequency) is greater than ZOO Hz. For this application, a prf up to 
10, 000 Hz can be used. A simple system to take advantage of this 
redundancy in the received pulse train is to integrate the received 
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Table 1 


Monostatic Radar Parameters 
3 cm 


\ 

10 cm 


30 cm 


100 cm 


size, dia. 
Antenna gain 

beamwidth 

Scattering volume 


4. 05 m 

13.5m 

15 m 

15m 

50 db 

50 db 

41. 4 db 

31 db 

. 519° 

. 519° 

1.4° 

4. 66° 

3. 86xl0' 2 

3. 86x1 0 

. 282 

3. 14 


R > 10 m R 

. 2 . -16 . -2/3 Q , 1n -15 

for C =10 cm 1.83x10 

n i 


R 


1. 83x10 


-15 


R R 

-15 -15 

1.27x10 .85x10 


NF 

RF amp. 

gain 

4 db 
15 db 

4 db 
15 db 

3 db 
15 db 

2 db 
15 db 

crystal mixer 
plus NF 
I. F. amp. 

9 db 

9 db 

9 db 

9 db 

Loss before 
R. F. amplifier 

. 5 db 

. 5 db 

. 5 db 

. 5 db 

T 

e 

760°K 

760°K 

500°K 

355°K 

T crystal mixer 

e • i 

receiver only 
(no RF amp) 

2350°K 

2350°K 

2350°K 

2 350°K 

Max. Range for 
detection of 

5. 5 Km 

18.5 Km 

21.2 Km 

20. 8 Km 


2 -16 -2/3 

C =10 cm 
n 


with P =10 watts 

s 

(single pulse— =1 
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signal with a low pass filter after the detector. For a 10 kHz prf and a 100 
Hz low pass filter, the number of pulses integrated is — = 100. The 
noise is uncorrelated from one pulse to the next and thus does not add in 
the integrator as the signal does. If the signal is detected in an envelope 
detector, the phase information is lost and the integration efficiency is much 
less than expected. Figure 19 shows the relationship between the number of 
pulses integrated and the integration improvement factor defined as 

— (after integration) = I— (single pulse)' . 

For a multiplier type detector as used with a coherent radar, 

I = n 

where n = number of pulses integrated. 

In the example above the 100 pulses integrated, 

S S 

effective— = 35— (single pulse) 

N N 

when an envelope detector is used. 

The transmitted pulse may be used to phase synchronize an 
oscillator which will then supply a reference signal of the proper phase to 
a multiplier type detector. The result will be that the phase information 
will be retained and the integration improvement factor will be equal to 
the number of pulses integrated, or 100 for the example. This may be 
done as shown in the block diagram in Figure 20. 





Block Diagram of a Coherent Pulse Radar 
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B. Bistatic Radar 

The bistatic radar performance when using a pulse system is 
substantially the same as the vertically pointing backscatter system except 
for the increase in a due to the reduced scattering angle. This increase 
will be about 5 db for a 90° scattering angle. However, the bistatic radar 
allows the scattering volume to be defined entirely by the intersection of 
the beams of the two antennas thus allowing CW operation. 

CW operation of a radar for the 100 cm wavelength for a trans- 
mitter power of about 1000 watts is relatively easy to obtain for low cost 
compared with a microwave system of the same power. Using this coherent 
transmitter ,1 the receiver can use a reference from the transmitter to 
mix with the received signal from the disturbed region. A multiplier 
type detector would probably be the best to use as it has a linear character- 
istic. By this means both the amplitude spectrum and doppler of the 
turbulent air can be measured and knowledge of the speed of the moving 
air can be obtained. 

j Moderate turbulence has a gust velocity of 20 to 30 mph or 9 to 
13 meters/sec. For vertical motion, the doppler frequency shift is 
given by 



V = vertical velocity of target* 


C = phase velocity. 
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Therefore, the doppler frequency for the 100 cm radar will vary between 
18 to 26 Hz. This means that a very small receiver bandwidth of about 
50 Hz can be used. 

By putting a 50 cps filter at the output of a linear detector, the 
effect of a 100 cps I. F. bandwidth can be obtained. Figure 21 illustrates 
the system for this case. 

With a coherent system of this type, two methods of clutter 
reduction are possible. The transmitter signal can be directly carried 
to the receiver by a coax cable and then adjusted in amplitude and phase 
to cancel the clutter; also, a separate antenna can be used to pick up the 
transmitter signal directly and this signal added to the signal from the 
main antenna to cancel the clutter. 

To obtain a comparison with the monostatic system, a repre- 
sentative maximum range calculation is made using the same parameters 
as before except for P^_ = 1 Kw and receiver bandwidth = 100 Hz. The 
antennas are assumed to be 2 Km apart, elevated 45° to give 9 = 90° 
and use horizontal polarization to give g = 90°. 

backscatter bistatic CW 

single pulse coherent detector 

4 / P A 4. 3 x 10 -4 . 086 . 872 

N t 

For a 10 db signal-to-noise ratio, the transmitted power would 
have to be about 25 Kw (pulse) for backscatter or 11.5 watts for the bistatic 



Block Diagram of cw Bistatic Radar 
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system (CW). The average power of the coherent pulse radar for a prf of 
10 kHz and pulse length of 2p,s would be 23 watts or twice that of the CW 
radar. This is due to the smaller scattering volume for the backs catter 
pulse radar. A precedent for this low power is the system used at Boeing 
using a 74 mHz, 100 watt CW transmitter and a receiver bandwidth of about 
50 Hz. Signals from disturbed refractive index at altitudes of 20, 000 to 
30, 000 feet of up to 10 db over the noise level. 

The main disadvantage is that separate antennas are required; 
however, this probably would be necessary for the backscatter radar in 
order to achieve the low noise and clutter levels required. 

C. Comparison with Other Systems 

Other researchers have detected radar echoes vhich are due to 
some type of disturbed refractive index of the atmosphere. Among these 
are Saxton, et al. Atlas, Hardy and Glover^; and Buehler and 

Lunden^’ ^ Atlas, et al, have detected turbulence associated 

with a tropopause at about 10 to 12 Km. The measured radar cross 

-17-1 2 -16 
section was about 7.7x10 cm corresponding to = 4. 4 x 10 

-2/3 

cm . Saxton, using a much less powerful radar, detected refractive 
index variations due to turbulence in layers at altitudes of 2 to 3 Km. 

These observations by Saxton were not measured quantitatively. 

Buehler and Lunden at a site in Seattle, Washington, have seen 
many times, with a VHF radar, echoes which are apparently associated 
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with turbulence. These echoes have been observed up to about 40, 000 feet 
and correlated with reports of turbulence in the area. ^ Another system 
was installed in a Boeing 727 jet and flown in areas of turbulence. A total of 
eight encounters with turbulence were recorded. The radar showed the 

(g\ 

presence of turbulence in each case' 

Table 2 gives the parameters of the Boeing system, and Wallops 
Island (Atlas, et al). 

Comparison of Table 2 with Table 1 indicates that the designed 

systems are not as good as the Wallops Island radar but better than the 

Bbeing systems. The Wallops Island radar has detected a return from 

near the tropopause at arrange of 12 Km. From these measurements, they 
2 -16 -2/3 ^ 

calculated C = 4.4 x 10 cm , which lends support to the 
n 

assumptions upon which the present system design is based. The Boeing 

ground system has detected disturbed refractive index at altitudes up to 

( 12 ) 

38, 000 feet with signal-to-noise ratios of 6 db or more. The Boeing 

airborne radar has detected turbulence and encountered it with the airplane 
eight times during April and May 1966. The detection range is equal to the 
altitude of the airplane due to ground return, but signals from light turbulence 

(g) 

(. 3g vertical acceleration) were observed at ranges up to 30, 000 feet 

Calculations and measurements with an M33 radar as reported in 
reference 23 show that the X-band radar will not detect turbulence. However, 
a few returns have been received that have not come from small targets but 


possibly from extended layers. 
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The characteristics of this radar are 

P T = 250 Kw 

r = . 25ns 

G = G = 40 db 
T R 

Noise figure = 13 db. 


2;lw 
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Table 2 

Existing Systems 

Boeing System Wallops Island 

ground air 

Transmitter 


X 

1. 38 m 

1.38m 

70 cm 

10.7 cm 

3. 2 cm 

P T 

120 Kw 

200 Kw 

7. 5 Mw 

4. 5 Mw 

170 Kw 

PRF (pps) 

60 to 200 

60 

960 or 320 

320 

320 

T 

10 (is 

2-6 |j.s 

. 1 to 6 . 2 |j,s 

2.2 |j.s 

1 . 8 (is 

Pav 


50 w 






Receiver 




NF 

4 db 

3. 5 db 


• 

3. 5 db 

T 

e 


420°K 

550°K 

500°K 


Bw 

200 kHz 


. 5 to 20 mHz 

650 kHz 

7 50 kHz 

Antenna gain 

22 db 

18. 3 db 

35 db 

51 db 

56 db 

Max. range for 

r 2 -16 - 2/3 

C =10 cm ° 

i ■ • - 

6 . 6 Km 

3. 9 Km 

61.8 Km 

75. 7 Km 

12. 76 Km 



VIII. CONCLUSIONS 


This study of the radar parameters indicates that disturbed refractive 
index in the atmosphere may be detected by means of radar measurements 
with moderate power at short ranges. The most important features of the 
radar are antenna size and a receiver sensitivity. It is concluded that there 
is little difference in effectiveness of the system for wavelengths from 1 m 
to 10 cm as long as the same antenna size is used and the antenna beam is 
filled. For the purpose of analyzing the disturbed refractive index to 
determine what type of return signifies danger to airplanes, the smaller 
scattering volume obtained with wavelengths of 10 cm to 30 cm would 
probably give a more accurate estimate of the structure of the turbulent 
area as it passes through the antenna beam. 

There are several practical difficulties which limit the usefulness of 
most present radar systems for turbulence detection. The most serious 
of these is minimum range performance. Recovery time and clutter level 
usually restrict most radar systems to a relatively long range at which the 
turbulent scattering signal is lost in the noise. 

Therefore, even though little success has been attained with conventional 
radar in detecting turbulent regions in the atmosphere, a specially designed 
research system should have little trouble in the detection and analysis of 
the signal from regions of disturbed refractive index of the atmosphere at 
lower elevation. 
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The optimum wavelength will lie sromewhere around 10 cm to 30 cm. 

This choice is due to the wavelength dependence of the scattering cross 
section and the limits on antenna size and gain. 

The antenna will have a maximum diameter of 15 meters and a maximum 
gain of 50 db. At wavelengths shorter than 11 cm, the size is limited by 
the gain consideration. 

In order to use moderate transmitter power (about one to three mega- 
watts peak), the receiver will need to have a low noise preamplifier. From 
these studies, it is concluded that a radar for airborne detection of disturbed 
refractive index has severe practical limitations. The most important 
limitations are antenna size and transmitter power. 
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